A real space decoupling transformation for quantum many-body systems 
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We propose a real space renormalization group method to explicitly decouple into independent 
components a many-body system that, as in the phenomenon of spin-charge separation, exhibits 
separation of degrees of freedom at low energies. Our approach produces a branching holographic 
description of such systems that opens the path to the efficient simulation of the most entangled 
phases of quantum matter, such as those whose ground state violates a boundary law for entangle- 
ment entropy. As in the coarse-graining transformation of [Phys. Rev. Lett. 99, 220405 (2007)], the 
key ingredient of this decoupling transformation is the concept of entanglement renormalization, or 
removal of short-range entanglement. We demonstrate the feasibility of the approach, both analyti- 
cally and numerically, by decoupling in real space the ground state of a critical quantum spin chain 
into two. Generalized notions of RG flow and of scale invariance are also put forward. 

PACS numbers: 03.67.-a, 03.65.Ud, 02.70.-c, 05.30.Fk 
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At low energies, an extended quantum many-body sys- 
tem may decompose into several independent compo- 
nents, as in the phenomenon of spin-charge separation, 
whereby a system of interacting electrons decouples into 
spinons and holons [1]. Another example is a Fermi liquid 
in D > 2 dimensions, whose Fermi surface factor izes into 
an infinite collection of one-dimensional critical systems 
Q. More exotic systems, such as spin Bose- metals dis- 
playing a Bose surface [3], are also believed to decouple 
in a similar way. In the above examples, the factorization 
of the many-body Hilbert space into components at low 
energies becomes evident when the Hamiltonian is ex- 
pressed in terms of the adequate collective quasi-particle 
operators in momentum space. However, a generic inter- 
acting many-body system is naturally modeled by some 
microscopic degrees of freedom with a local Hamiltonian 
in real space, where this factorization may no longer be 
manifest. 

The purpose of this paper is to propose a theoreti- 
cal framework and a numerical technique to obtain, in 
real space, a factorized description of the ground state 
(more generally, low energy subspace) of a many-body 
system that experiences separation of degrees of freedom 
at low energies. There are several good reasons, both 
conceptual and computational, to search for one such de- 
scription. Fermi liquids and spin-Bose metals, which are 
among the most entangled phases of quantum matter, 
pose an exceptional challenge to real space renormaliza- 
tion group (RG) methods. Recall that in D > 2 di- 
mensions, where the ground states of many phases (both 
gapped and gapless) conform to a boundary law for en- 
tanglement entropy [4] , the ground states of systems with 
emergent Fermi or Bose surfaces display logarithmic vio- 
lations of this boundary law [5] . A logarithmic correction 
may seem small business. However, it has the dramatic 
effect, as one flows towards large length scales, of increas- 
ing the number of degrees of freedom that one needs to 



deal with, thus rendering real space RG methods inef- 
ficient. In particular, it precludes the use of all known 
forms of tensor network ansatze Q, which can only 
efficiently describe ground states that strictly obey an 
entropic boundary law. This is most unfortunate, since 
tensor networks have recently emerged both as the ba- 
sis of powerful non-pertubative approaches to many-body 
systems on the lattice and as the natural language 
to describe and classify exotic phases of quantum matter 
[Tq| . An ultimate goal of our work is to understand the 
entanglement pattern in systems, such as Fermi Liquids 
and spin-Bose metals, that violate the entropic boundary 
law, and to extend the tensor network formalism to these 
exceptionally entangled phases of quantum matter. 

Specifically, here we will show that entanglement 
renormalization, or removal of short-range entanglement 
— already the basis of a renormalization group transfor- 
mation for quantum systems on a lattice — , can also 
be used to identify and decouple, in real space, differ- 
ent components of a system that experiences separation 
of degrees of freedom. For illustrative purposes, we will 
demonstrate the performance of the approach in D = 1 
dimensions only, by numerically decoupling a critical 
quantum spin chain into two independent components. 
However, most of our discussion also readily applies to 
D > 2 dimensions. Decoupling in real space leads to 
a generalized multi-scale entanglement renormalization 
ansatz (MERA), the branching MERA [Hj, which is a 
tensor network ansatz capable of reproducing violations 
of the entropic boundary law. It also defines a holo- 
graphic geometry [13] with branches for each of the com- 
ponents, drastically reducing computational costs, and 
suggesting a generalized notion of scale invariance. 

Let £ be a one-dimensional lattice made of N sites, 
where each site is described by a Hilbert space V of finite 
dimension x- We denote by H and |\&) a local Hamil- 
tonian acting on Y® N and its ground state, respectively. 
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FIG. 1. (color online) Decoupling transformation V for a 
one- dimensional lattice C. V maps the lattice model into 
two decoupled lattice models, with lattices Ca and Cb- The 
decoupling transformation V decomposes as the product of 
two rows of (two-site) disentanglers u and v, and a number of 
(two-site) swap gates that bring all the sites of sublattice Ca 
(Cb) together to the left (respectively, right). 
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Here we are interested in a very particular class of lattice 
models. Namely, we will assume that the model can be 
decoupled into two models, with lattices Ca and Cb, each 
made of N/2 sites, under a suitable decoupling transfor- 
mation V to be introduced below. 

Let us first clarify what we mean by 'decoupling' a 
lattice model into two, since this term may be used to 
describe several inequivalent scenarios. The strongest 
notion refers to 'Hamiltonian decoupling': V maps the 
Hamiltonian H into a sum of two Hamiltonians Ha and 
Hb for sublattices Ca and Cb, 

H -+H' = V ] HV = H A + H B - (1) 

Notice that the ground state is then mapped into the 
product of ground states for Ha and Hb, 

\*)^\*') = V*\*) = \*a)®\*b)- (2) 

A second, weaker notion of decoupling is that of 'ground 
state decoupling', where the transformed Hamiltonian H' 
may contain interactions between sublattices but, never- 
theless, its ground state \^') still factorizes as in Eq. O 
We note for later reference that in this case \S&a) and 
\^b) are also the ground states of the reduced Hamilto- 
nians H' A and H' B , 

H' A = {^ B \H'\^ B ), H' B = {^ A \H'\^ A ). (3) 

Decoupling transformation. — Our decoupling transfor- 
mation V is a unitary transformation implemented in real 
space by means of local gates, as depicted in Fig. [TJ It is 
made of two types of local gates: (i) disentanglers, that 
in the present example appear organized in two rows; 
and (ii) swaps, used to regroup the sites of sublattices 
Ca and Cb- Disentanglers u and v are two-site transfor- 
mations, i.e. u : V® 2 V® 2 and v : V 02 V 02 , con- 
strained to be unitary, i.e. uv) = v)u = vv^ = v^v = I 02 , 



FIG. 2. (color online) (a) Decoupling transformation V, 
mapping the ground state of a local Hamiltonian H on 
lattice £, into the tensor product of states \^a) and \^b) 
for lattices Ca and Cb, as in Eq. [2] (b) Coarse-graining 
transformation W, mapping the ground state |^) of a local 
Hamiltonian H on lattice C into the state \^') of a coarse- 
grained lattice C! . 



with 0(x 4 ) complex coefficients, which can be treated as 
variational parameters (to be determined, for instance, 
through numerically optimization). The role of disentan- 
glers is to remove short-range entanglement between sub- 
lattices Ca and Cb so as to produce a factorized ground 
state, Eq. [2] and Fig. [2fa), or to eliminate interaction 
terms that couple the two sublattices to produce a decou- 
pled Hamiltonian, Eq. [TJ Disentanglers are then followed 
by swap gates, each of which permutes the basis states 
of a pair of lattice sites. For bosonic degrees of freedom 
a swap gate acts just as 

swap : \i) (g) \j) — » \j) g) (bosons) (4) 

where {\i}} denotes a basis in V, whereas for fermionic 
degrees of freedom it includes a conditional minus sign, 

swap : 

K) ® |j) -> \j) <8> \i), (fermions) (5) 

where f(i,j) = 1 if both \i) and \ j) carry an odd fermionic 
particle number and f(i,j) = otherwise. For anyonic 
degrees of freedom, the swap gates will embody the ap- 
propriate representation of the braid group. 

A key question is what type of many-body systems are 
susceptible to decouple, in the sense of Eqs. [T]or[2j under 
the decoupling transformation V. A concrete example is 
given by the XX quantum spin chain, 

F xx = J2 (Mr+l + Y r Y r+1 + fi/2Z r ) , (6) 
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FIG. 3. (color online) (a) Sequence of coarse- graining trans- 
formations W and decoupling transformation V for a system 
with separation of degrees of freedom, with s* = 3. Each 
transformation V or W decomposes into two-site gates as in 
Fig. [2] The resulting quantum circuit defines the branching 
MERA. (b) Corresponding holographic tree. The decoupling 
transformation V produces branching in the holographic or 
scale direction (vertical axis). 



where X, Y and Z are Pauli matrices. As discussed in 
Ref. [14] (exercise 12.1, pp 480), for \i = this Hamil- 
tonian can be exactly decoupled into two independent 
copies of the quantum Ising chain with critical transverse 
magnetic field, 



Rising = (X r X r+1 + Z r 



(7) 



Moreover, as can be easily seen by rewriting i^ xx in Ma- 
jorama fermionic modes jl5[, this Hamiltonian can be 
decoupled as in Eq. [I] by a transformation V made of an- 
alytical disentanglers (that simply permute some of these 
modes) and the fermionic swaps of Eq. [5] (see Appendix). 

Holographic branching. — Generically, however, sepa- 
ration of degrees of freedom is likely to take place at a 
length scale, the decoupling length A*, that is much larger 
than the microscopic lattice spacing, i.e. A* ^> 1 [we ex- 
press length scales in units of the lattice spacing in C). 
In this case, we will proceed by first coarse-graining the 
microscopic model into an effective lattice model whose 
lattice spacing is of the order of A*. 

Let W denote the coarse-graining transformation of 
Ref. [11], shown in Fig. Efb), which maps a lattice C 
with N sites into an effective lattice £ with N/2 sites, 
and which uses disentanglers u in order to remove short- 
range entanglement, before applying isometries w to map 
two sites into one. [Those familiar with Refs. may 
have noticed that W is a particular case of V in which 
one of the two resulting lattices is in a product state]. 
We will use the scale parameter s = log 2 (A) to label the 
lattice obtained after applying s times the trans- 
formation W on C. Then, in a system with decoupling 
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FIG. 4. (Colour online) (a) Scaling dimensions for scaling 
operators (extracted from a scale- invariant MERA [8|]) of each 
of the two decoupled spin chains. They are organized accord- 
ing to the conformal towers [14] of the critical quantum Ising 
model, (b) The entanglement entropy in each of the two de- 
coupled chains (central charge c — 1/2) is half of that in the 
system without decoupling (central charge c = 1), see Eq. [8] 



scale s* = log 2 (A*), we obtain a factorized description 
by: (i) applying a total of s* coarse-graining transforma- 
tions W to the microscopic model, which result in a se- 
quence of effective lattices {C^°\C^\ • • • , C^ s )}; (ii) ap- 
plying the decoupling transformation V on C^ s \ which 

and 



produces two decoupled models with lattices C 
£g and (iii) applying transformations W individu- 
ally on each of the decoupled models to further coarse- 
grain them, as illustrated in Fig. [3l 

Recall that the coarse-graining transformation W 
alone produces the MERA [3, [ll[ . In much the same way, 
if we collect all the two-site gates (disentanglers, isome- 
tries and swaps) contained in transformations W and V 
in Fig. [3j we obtain a new class of tensor network state, 
which branches into two independent tensor networks at 
scales s larger than the decoupling scale and which, 
accordingly, we refer to as the branching MERA [T3 |. 
This is a variational ansatz for states of the microscopic 
lattice £(°\ which inherits many of the computational 
advantages of the MERA, such as the efficient evaluation 
of local expectation values, e.g. energy and local order 
parameters, and of correlation functions. 

To illustrate the performance of the approach, let us 
return to the critical quantum spin chain of Eq. [6] For 
j(i / 0, Hamiltonian i^ xx can no longer be decoupled 
analytically into two copies of Hamiltonian H Isins , be- 
cause the magnetic field ^ Z r (or chemical potential 
in fermionic variables) turns into an interaction between 
the two sublattices. However, we managed to numerically 
factorize the ground state of the model. For instance, 
for \± = 2cos(7r/4) (or 1/8 filling in the fermionic rep- 
resentation), an infinite lattice was first coarse-grained 
three times (i.e. s* = 3) with transformation W and 
then decoupled into two chains with a transformation V 
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FIG. 5. (color online) (a) Holographic tree without branch- 
ing, representing the coarse-graining of a critical lattice model 
without separation of degrees of freedom, (b)-(c) Other pos- 
sible holographic trees, with a regular pattern of branching. 
They describe fixed-points of a (generalized) RG flow where 
separation of degrees of freedom occurs at all length scales. 



containing fermionic swaps. Then, assuming a factorized 
ground state (Eq. [2]) at that length scale, we used a scale- 
invariant MERA to find the ground state of each of the 
two reduced Hamiltonians (Eq. [3]). All variational pa- 
rameters of the branching MERA were numerically op- 
timized with simple variations of the standard, energy 
minimization algorithms of Refs. 0, [l6| . 

Specifically, using a local dimension ^ = 16 (x = 8) 
for scales s smaller than (respectively, larger than) the 
decoupling scale s* = 3, the energy optimization took 
less than a day on a 2.66GHz workstation with 8GB of 
RAM. The resulting branching MERA approximated the 
exact ground state energy per site with an accuracy of 
8 x 10 -7 and duly captured the structure of large scale 
correlations in the system: from each branch we were able 
to extract [8[ accurate conformal data corresponding to 
the conformal field theory [14] that describes the critical 
Ising chain of Eq. [7J see Fig. 2J 

Discussion. — Two remarks are in order. First, in 
the example above decoupling was possible by using 
fermionic swaps (any other choice of particle statistics re- 
sults in persistent interactions/correlations between the 
two sublattices), and applying V at a proper scale s*. 
The particle statistics and decoupling scale s* constitute 
important knowledge which we acquired by trial and er- 
ror. Second, let pL denote the ground state reduced 
density matrix for a block of L contiguous sites of £(°). 
The entanglement entropy Sl = — tr[pL log 2 (pz,)], which 
scales as Sl ~ | log(L) [15j, receives contributions from 
all length scales A < L. Consider a block L larger than 
the decoupling length A*. Then it can be seen that the 
contribution to Sl from a length scale A G [A*, L] is split- 
ted evenly in contributions from the two decoupled chains 
(each of which asymptotically contributing | log(L), Fig. 
[4]) which implies 

S L ^S L A) +S L B \ L«\* (8) 

In other words, the (real-space) entanglement of a large 
block is roughly the sum of entanglements from each of 
the components into which the model decomposes. 



Generalizations to D > 2 dimensions are best dis- 
cussed by introducing a holographic tree, Fig. [3fb), which 
captures the branching structure of a system under real 
space coarse-graining. In a general setting, we can imag- 
ine a holographic tree with a complex branching struc- 
ture, corresponding to several instances of factorization 
of degrees of freedom occurring at various decoupling 
scales. This leads to an intriguing, multi-system gen- 
eralization of the RG flow, with a variable number of in- 
dependent many-body systems as a function of the scale 
parameter s. During the flow towards low energies, new 
systems may be added through branching. In particular, 
a new notion of RG fixed point arises, corresponding to a 
self-similar holographic tree, see Fig. [5j which represents 
a many-body system with branching at all length scales. 

A branching MERA in D = 2 dimensions can describe 
the ground state of free fermions on a square lattice with 
a Fermi surface [l2j ]. The system seems to correspond 
to one of these novel, multi-branched RG fixed points. 
The entanglement entropy of a large block of L x L sites 
receives contributions from all relevant branches, general- 
izing Eq. [8] in such a logarithmic violation of the bound- 
ary law, S l ~ Llog 2 (L), is reproduced. Accordingly, we 
expect the branching MERA to also provide an efficient 
description of Fermi liquids and spin-Bose metals. 
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APPENDIX 

In this appendix we show that the XX quantum spin 
chain can be exactly decoupled into two (critical) quan- 
tum Ising spin chains by applying local disentanglers - 
provided that fermionic degrees of freedom are used to 
separate the two systems. This derivation follows from 
the answer to exercise 12.1 in page 480 of Ref. [lij . 

The XX quantum spin chain Hamiltonian without 
magnetic field is 

H = Y J XrX r+1 +Y r Y r+1 (9) 

r 

Let us introduce Majorana fermion operators c r and d r , 
given by 

C r = [l[Zl)Xr, d r = (n^) y - ( 10 ) 
\Kr J \l<r J 

Notice that all c r 's and <i r 's anticommute pairwise, ex- 
cept with themselves (since they fulfill (c r ) 2 = (d r ) 2 = 
I). In terms of these Majorana operators the Hamilto- 
nian reads 

H = i (c r d r +i — d r c r +i) . (11) 

r 

Further, for each odd r, we apply a unitary gate u on the 
pair of sites (r, r + 1), such that 

u c r = c r+ i, u d r v) = c r , 
u c r+ i v) = d r , u d r +\ v) = — d r +i. (12) 

[Notice that u preserves the anticommutation relations 
and it is therefore a canonical transformation of the 
fermionic variables.] Then the transformed Hamiltonian 
reads 

H = y ^ /i r?r +2? ^r,r+2 = * (<i r C r +2 — C r d r ) , (13) 
r 

or H = Ha + Hb, where Hamiltonians Ha and Hb, 

Ha = hr,r+2, H B = hr,r+2 (14) 

odd r even V 

commute with each other, since they are made of terms 
that are quadratic in Majorana fermion operators and 
they act on two different sets of sites. Our decoupling 
transformation V in Fig. [T]will then consist of a row of 
disentanglers u as in Eq. [12j followed by trivial disen- 
tanglers v (i.e. v = I) since the system has already been 



decoupled. To group together all sites of sublattice Ca 
to the left, and all sites of sublattice £^ to the right, 

(a) 




FIG. 6. (Colour online) (a) Representation of the exact dis- 
entangler u of Eq. [12] used to decouple Hamiltonian H in Eq. 
[9] Each site, represented by an oval, contains two Majorana 
fermions, depicted as brown boxes. The disent angler u per- 
mutes the Majorana modes and adds a minus sign to one of 
them, (b) Hamiltonian terms before (below) and after (above) 
the disent angler u is applied on pairs of sites. Each Hamilto- 
nian term, e.g. cicfe, is depicted with a line connecting two 
Majorana modes. Notice that after applying the disent an- 
gler, which e.g. maps c\d^ to — C2C/2, the two sublattices are 
no longer connected by interactions. 

we will use fermionic swap gates, Eq. [5] Notice that 
these swaps are trivial in fermionic variables, simply in- 
terchanging (c r ,d r ) with (c r+ i, d r +i). However, they are 
non-trivial when expressed in spin variables as in Eq. 
Once the two lattices have been separated, we can de- 
fine independent Jordan Wigner transformations on each 
lattice & a \ a = A,B, We can now define odd and even 
spin variables by 

4^ = (U Z i a) ) Xr"^ ^ = (U Z i a) ) Y r a) ' ( 15 ) 
\Kr J \l<r J 

in terms of which we obtain Hamiltonians 

= (X^xl% + ZW) , a = A,B, (16) 

r 

which correspond to two copies of the quantum Ising 
model with critical transverse magnetic field. 



